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1 Introduction 

The social mixing structure of a population or a group of interacting populations 
play a crucial role in the dynamics of a disease transmission || . The most well-known 
examples of epidemics which are spread by means of the interaction between different 
populations are those related to venereal diseases. These diseases are transmitted 
by sexual contacts between two or more different populations. An epidemiological 
model which considers these interactions is called multigroup model. A large number 
of multigroup models have been described in J7). These models usually lead to 
some high-dimensional systems of differential equations with a probably complicated 
dynamics |TJ| . In order to avoid these complications, we consider a single group and 



we assume that there is a disease transmission from the outside. This assumption 
is natural, since we are living in an open environment with a limited control on 
the outside. The effect of this assumption has been studied by using numerical and 
statistical methods. For example in 0, it has been claimed that most of Swedish 
HIV-patients with heterosexual transmission had been in contact with a partner of 
foreign extraction. The effect of the immigrant subpopulation has been considered 
by many others too |], |9], pLlfl . A statistical approach to this concept in |10j showed 
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that the subpopulation with a high prevalence of HIV-infected individuals is the 
immigrant one. 

In this paper, we consider a simple SIRS type model of disease transmission 
in an open environment. In an SIRS model, we divide the population into three 
classes consisting of Susceptible, Infective and Removed individuals. By an open 
environment we mean that our population has contacts with other populations and 



there is a disease transmission from the outside ||15|| . We assume that a proportion 
of susceptible individuals are infected in this way. We also assume that there is 
a special control on the people who are frequently in contact with foreigners and 
a proportion of those who are infected in this way are known and removed. In 
our model, the incidence function is of proportionate mixing type introduced by 
Nold [IS|. The demographic assumptions are also very simple. Natural births and 
deaths are assumed to be proportional to the class numbers with all newborns 
susceptibles and the excess deaths due to the disease among infectives and removeds 
are proportional too. We could have considered more complicated demographic or 
vertical transmission assumptions 0, f| , but these parameters have no mathematical 
significance and one can easily conclude that our main results would still be valid. 
We want to avoid these complications for two reasons. The first one is that even in 
this simple case, multiple endemic equilibria may occur. Indeed we shall show that 
for some suitable values of parameters, the proportions system admits two sinks and 
one saddle point in the feasibility region. The existence of multiple equilibria or limit 



cycles has been shown for more complicated systems such as multigroup models [14 



or a single group with nonlinear incidence functions [12, 16, 17 1. The second reason 
is related to the technique used here to determine the number of endemic equilibria. 
We state a simple model in order that we can exhibit our technique more clearly. 
The technique used here is based on a careful choice of Jordan curves and counting 
the number of rest points inside them. This technique has no hard analysis and can 
be easily applied to other similar systems |2~1| , [22| . The reader can verify that our 
results hold for similar SIR and SIRI systems as well. 

We first in the next section, state the model and some results concerning the 
non-existence of certain types of solutions. In Section 3, we use Poincare index to 
obtain some partial results in a special case. Then in Section 4, we use these results 
together with a perturbation method to examine the general case. Our analysis 
relies on an index lemma concerning the Poincare index of a class of Jordan curves. 
We provide the proof of this result at the end of this paper. 
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2 The Model 



We consider a model of disease transmission in as nonconstant population of size N 
divided into three classes: susceptibles, infectives and removeds, the number of each 
class is given by S, I, R, respectively. We set N = S + I + R and use the following 
parameters which are assumed to be positive unless otherwise specified: 

b= per capita birth rate, 

d= per capita disease free death rate, 

Si— excess per capita death rate of infectives, 

62= excess per capita death rate of removeds, 

A= effective per capita contact rate of infectives, 

a= per capita removal rate of infectives, 

7= per capita recovery rate of removeds. 

Here as mentioned before, we assume that there is a disease transmission from the 
outside to our population. Let (3 > be the per capita transfer rate of the disease 
into this population. Furthermore the susceptible individuals which are infected in 
this way, enter to the classes / or R of proportions (3\ and fa where fa + fa = fa 
The above hypotheses lead to the following system of differential equations in R^_, 
where " ' " denotes the derivative with respect to t, the time. 

( XI 9 

S' = bN-(d + P)S + jR-—, (2-1) 

V = ftS _ (d + e x + a)I + (2-2) 



R' = {3 2 S-(d + e2 + i)R + aI, (2-3) 



XI 



where — is of the proportionate mixing type introduced by Nold ||18|| . The total 



population equation is obtained by adding the above three equations: 

N' = (b - d)N - e x I - e 2 R. (2-4) 

SI R 
Now if we set s = — , i = — and r = — , the equations (2-l)-(2-3) yield 

= 6- (6 + /9)s + 7r + (^ - X)is + e 2 rs, (2-1)' 
= flxs - (b + ex + at)i + Xis + e x i 2 + e 2 ir, (2-2)' 
= fas- (b + s 2 + l)r + ai + Exir + E 2 r 2 . (2-3)' 
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We will determine the asymptotic behaviour of the solutions of this system. Prior 
to this we have the following concepts of ODE's related to our system. 

Given an autonomous system of ordinary differential equations in M. n 

d i = f{x) (2 - 5) 

We will denote by x.t the value of the solution of this system at time t that is x 
initially. For V C W\ J C R, we let V.J = {x.t : x e V,t e J}. The set V is 
called invariant if VM = V and it is called positively invariant if VM + = V. For 
Fci" the w-limit set of Y is defined to be the maximal invariant set in the closure 
of y.[0, oo). A closed curve connecting several rest points whose segments between 
successive rest points are heteroclinic orbits is called phase polygon. By a sink we 
mean a rest point at which all the eigenvalues of the linearized system have negative 
real parts. A rest point is called a source if these eigenvalues have positive real parts 
and it is called a saddle point if some of these eigenvalues have positive real parts 
and the others have negative real parts. A rest point is called nondegenerate if all 
of these eigenvalues are nonzero and it is called hyperbolic if all of its eigenvalues 
have nonzero real parts. 

Now we continue the analysis of the system (2—1)'— (2— 3)'. If we set = s+i+r 
then ^2' = (1 — ^2)(b — E\i — Sir). Therefore the plane = 1 is invariant. We 
consider the feasibility region 

D = {{s,i,r)\s + i + r = 1, s > 0, % > 0, r > 0}, 

which is a triangle and on its sides we have 

s = 0^s' = b + 'yr , % = =^ i' = fas , r = r' = f3 2 s + ji. 

Since all parameters are positive, D is positively invariant and any solutions of the 

o 

system (2 — 1)' — (2 — 3)' with initial point in dD immediately enters D, where 

o 

dD and D are the boundary and the interior of D, respectively. From now on, we 
examine the dynamics of this system in the feasibility region D. 

Using the relation s + i + r = 1, we see that the system is essentially two 
dimensional. Thus we can eliminate one of the variables to arrive at the following 
quadratic planner system: 

s' = b + ~f + (e 2 — b — (3 — 7)s — 72 — e 2 s 2 + (si — e 2 — X)is, (2 — 6) 
i> = p lS + ( e2 - £l - a - b)i + (A - e 2 )is + (d - e 2 )i 2 . (2 - 7) 
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Notice that this planner system has at most four rest points. Moreover the dynamics 
of the system (2 — 1)' — (2 — 3)' on D is equivalent to the dynamics of this system 
in the positively invariant region D\ = {(s, i) \s > 0, % > 0, s + % < 1}. 

The following theorem is a special case of the results of || concerning the non- 
existence of certain types of solutions. 



Theorem 2.1. Let in (2-5), / be a smooth vector field in M. 3 . Let T(t) be a closed 
piece-wise smooth curve, which is the boundary of an orientable smooth surface 
Scl 3 . Suppose g : U — > IR 3 is defined and smooth in a neighborhood U of S with 
g(T(t)).f(T(t)) > and curl(g).n < 0, where n is the unit normal to S. Then V is 
not a finite union of the orbits of the system (2.5). 



Proposition 2.2. The system (2 — 1)' — (2 — 3)' has no periodic orbits, homoclinic 

o 

orbits, or phase polygons in D. 

Proof. In order to apply the above theorem, we define g = g\ + g-i + g$ where 

h{h r ) h(hr)' 



9i[},r) 
92(s,r) 
gs(s,i) 



0. 



If 



0, 



ir 
fi(s,r) 



sr sr 
f2(s,i) fi(s,i) 



SI 



SI 



,0 



and fx, fi and f 3 are the right hand side of (2 — 1)', (2 — 2)' and (2 — 3)' re- 
duced to functions of two variables by using ^2 = 1 respectively. Now after some 
computations |12|, we get 

curhj(s,z,r).(l,l, 1) = - 



ri 2 ir 2 sr 2 



is* rs* 

Since here g.f = 0, the proof is complete by Theorem 2.1. □ 



Corollary 2.3. The w-limit set of any orbit of the system (2 — 1)' — (2 — 3)' with 
initial point in D is a rest point. 

Proof. Since the vector field related to the system (2 — 1)' — (2 — 3)' is inward 
on 3D and D is compact, the cu-limit set of each orbit with initial point in D, is 

o 

a nonempty subset of D. By generalized Poincare-Bendixon theorem |2(| and 
Proposition 2.2., this set must be a rest point. □ 
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Corollary 2.4. The system (2 — 1)' — (2 — 3)' has no source in D. 

o 

Proof. Suppose there is a source in D for this system. Since D is positively 
invariant and there are finitely many rest points in D (at most four, since our 
system is quadratic), there must be infinitely many heteroclinic orbit running from 

o 

this source to another rest point. So there is a 2-gons in D which is impossible by 
Proposition 2.2. □ 

o 

Corollary 2.5. Every nondegenerate rest point of the system (2-6), (2-7) in D 1 is 
hyperbolic. 

o 

Proof. Let L be the linearization of this system at a rest point (s*,i*) in D 1 . We 
have to show that if det L > 0, then trace 1^0. We compute trace L at a rest 
point. 

{ds' 
— = (e 2 - b - /3 - 7) - 2e 2 s + (e x - e 2 - X)i, 
di' 
-q7 = (£2 + £1 - ol - b) + (A - e 2 )s + 2(ei - e 2 )i. 

From i' — and s' — at (s*, z*), we have 

trace L = — e 2 s h (£1 — e 2 )i ■ 

s* i* 

If trace L = 0, then we can slightly increase S\ to get trace L > and determine 
A and a so that (s*,i*) remains a rest point for the new values of parameters. We 

o 

may also assume that det L > at this rest point. Thus we obtain a source in D\ 
which contradicts Corollary 2.4. □ 

Remark 2.6. A nondegenerate rest point of the system (2 — 6), (2 — 7) is obtained 
by a transversal intersection of two conic sections s' — and i' = 0. In Section 4, 
we shall prove that this intersection is almost always transversal. 



3 A Special Case 

In this section we consider the planner system (2-6), (2-7) in the case b — f3\ — 7 = 0. 
These assumptions yields the following system of equations 

s' = s(e 2 - (3 - e 2 s + (e l - e 2 - X)i), (3 - 1) 

i' = i(e 2 - Si - a + (A - e 2 )s + (e 1 - e 2 )i). (3 - 2) 
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First of all notice that there are two invariant lines s = and i — with three rest 



points (0,0), (0,1 - 
at (0, 0) is 



) and (1 ,0). The matrix of the linearized system 

£2 — £1 £2 



£2-/? 
£2 — £1 — a 



with the eigenvalues To := £2 — P and Ti := £2 — £\ — a. 



The matrix of the linearized system at (0, 1 
e 2 -(3 + (e 1 -e 2 -\)(l- 



a 



£2 - £1 



)is 



a 



(X-b-e 2 )(l 
with the eigenvalues — Ti and 



£2 - £1 



£2 - £1 



) 



£1 - £ 2 + a 



T 2 :=e 2 -(3 + 



Notice that the rest point (0, 1 — 



(£1 - £2 - A)(e 2 - £1 - a) 



a 



£2 - £1 

belongs to Di if and only if T\ > and 



£2 - £1 

coincides with (0, 0) in the case of equality. It is easy to see that when 7\ < 0, the 

(X 

origin attracts the segment D\ n {s = 0} and if T\ > 0, then (0, 1 ) attracts 



(i>in{ a = 0})-{(0,0)}. 



£2 - £1 



The matrix of linearization of (3 — 1), (3 — 2) at (1 , 0) is 

£2 



P-e 2 


with the eigenvalues — T and 

T 3 : 



( £l _ e2 _A)(l- A) 

£2 

(£ 2 -£i-a)-(£ 2 -A)(l-^-) 

£2 -I 



£2 — £1 — a — 



/3 



(£2 -A) (£2-/5) 



£2 



Notice that the rest point (1 , 0) belongs to D\ if and only if T > and 

£2 

coincides with (0, 0) in the case of equality. It is easy to see that when T < 0, the 
origin attracts the segment D\ D {i — 0} and if Ti > 0, then (1 — — , 0) attracts 

(An = 0}) -{(0,0)}. 
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Proposition 3.1. If T < and 7\ < 0, then the origin is the only rest point of 
the system (3-1), (3-2). 

a (3 

Proof. In this case (0, 0) is a sink and the rest points (0, 1 ) and (1 , 0) 

£2 — £i £2 

are outside of D\. Since our planar system is quadratic, there are at most one rest 

o 

point in D\. Moreover if such a rest point exists, it must be nondegenerate, hence 
its Poincare index is ±1. For a Jordan curve C, let //+(C) and /x_(C) denote the 
number of rest points inside of C with the Poincare index +1 and -1 respectively. 
Here let C be the curve as shown in Figure (3-1) which is the boundary of the union 
of Di and a small disk centered at the origin. Then on this Jordan curve, our vector 
field is always tangent or inward. It follows from Lemma 5.1. (cf. Section 5) that 
Ic(X) = 1 where X is the vector field related to the system (3-1), (3-2) and Ic(X) 
is the Poincare index of X with respect to C. Now we can use the Poincare theorem 
to obtain fi+(C) — fJ>-(C) = 1. Moreover we have shown that £t+(C) + ft -(C) < 2. 
Therefore fi+(C) = 1 and fi-(C) = which finishes the proof. □ 




Figure (3-1) 
The Jordan curve related to Proposition 3.1. 



Proposition 3.2. If T > 0, T x > 0, T 2 < and T 3 < 0, then the system (3-1), 

o 

(3-2) has a saddle point in D\. 

ol 

Proof. These assumptions mean that the origin is a source and (0, 1 ) 

£2 - £1 

and (1 ,0) are sinks. It is also concluded that the system (3-1), (3-2) has only 

£2 

nondegenerate rest points. Let -Bi,i? 2 and B3 be small disks centered at the above 
three rest points respectively. Let C = dA and A = (Di U B 2 U B 3 ) — B\. (See 
Figure (3-2).) The Poincare index of the Jordan curve C is 1 by Lemma 5.1. and 
we can use Poincare theorem to obtain /x+(C) — fJ,-(C) = 1. Moreover we have 
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/i + (C) + /x_(C) < 3 and /i+(C) > 2. Therefore /x + (C) = 2 and /z_(C) = 1 which 

o 

means that there is a saddle point in D\ . □ 




Figure (3-2) 
The Jordan curve related to Proposition 3.2. 



Remark 3.3. The assumptions of the above proposition do not contradict. To see 
this suppose that T > and T\ > 0. Then it is easy to check that 

T 2 <0^<1 + ^— , T 3 <0^>(l-A)-i 

J-l e 2 — £\ J-l ^2 

(l-±r 1 <(l + -±-)^e l >\. 

Now if e 2 > Si > A, one can choose a > and (3 > so that T > 0, T\ > and 
T A A 

— € ((1 1 H ) to satisfy the assumptions of Proposition 3.2. This 

T\ e 2 £2 — ^i 

would be helpful for the reader who is more interested in numerical simulations. 



4 The General Case 

In this section we investigate the dynamics of the proportions system (2 — 1)' — (2— 3)' 
in the general case. In order to do this, we discuss the existence and stability of the 
rest points of the planner system (2.6), (2.7). Recall that the feasibility region D\ 
is positively invariant and this system has no rest point on dD\. Indeed the vector 
field corresponding to the system (2.6), (2.7) is strictly inward on dD\. Thus by 
Lemma 5.1., the Poincare index of this vector field with respect to dD\ equals 1. 
This is the first step to prove the main result of this paper. 
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Theorem 4.1. If the system (2-6), (2-7) has only nondegenerate rest points in 

o 

Di, then one of the following statements holds: 

o 

(A) There exists a unique rest point Di which is a sink and attracts Di(GAS). 

o 

(B) There are two sinks and a saddle point in D\. 

Moreover both of them occur for suitable values of the involved parameters. 

o 

Proof. All rest points in D\ are hyperbolic by Corollary 2.5. Let [i\ and \ii be 

o 

the number of sinks, saddles and sources in D\. Since the index of dDi is 1, we have 
A*o ~~ A*i + A*2 = 1 va. Di. Furthermore fj, + Hi + /x 2 < 4 and // 2 = by Corollary 2.4. 
Thus we have either fi = 1, ji\ = /i 2 = or /i = 2,/ii = l,/i 2 = 0. Now the first 
conclusion gives (A) and the second one gives (B). In order to see that each of the 
two above statements occurs, recall that in the case of Theorem 3.1., there is neither 
a saddle nor a nonhyperbolic rest point in D ± . Thus for small values of b\ and 7, 

o 

there cannot be any saddle point or nonhyperbolic rest point in D\. This means 
that (A) occurs. Similarly, under small perturbation, the saddle point obtained in 

o 

Theorem 3.2 remains in Di and yields the case (B) of this theorem. □ 

Remark 4.2. It is well-known that the basin of attraction of a sink is a connected 
open set. In the case (B) of the above result, Di is the distinct union of the basins 
of attraction of these two sinks and the stable manifold of the saddle point in D ± . 
The basins of attraction are open and D 1 is connected. Thus the stable manifold of 
the saddle point separates them. In order to specify these sets, we can numerically 
find those two points at which a>limit set changes from one sink to another when 
someone moves on dDi. (See Figure 4.1.) 




Figure (4-1) 

The phase portrait of the case (B) of Theorem 4.1. 
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Remark 4.3. The non-uniqueness of endemic equilibrium proportions yields some 
interesting conclusions. The most significant one is that the initial condition of the 
population may also be important besides the involved parameters. The effect of 
the initial condition is more crucial when the population equations (2-l)-(2-4) is 
considered. From (2-4) we get 

N ' h A ■ 

— = b - d - eii - e 2 r. 

Now suppose that a solution (s(t), i(t), r(t)) of the system (2 — 1)' — (2 — 3)' tends to 

b 

an equilibrium (s*, i*, s*). If we set T = , then N(t) — > oo if T > 1 

d + £\i* + £ 2 r* 

and N(t) — > if T < 1. (See || for more details.) Now each endemic equilibrium 
gives a T and when there are two endemic equilibria, we may get different values 
for T at these two points. 

It remains to consider the case in which there is a degenerate rest point in Di for 
our planar system. Let Q be the parameters space of the system (2-6), (2-7) as an 
open subset of IR+ and Qi be the set of all possible values of parameters for which 
the system (2-6), (2-7) has a nonhyperbolic (or equivalently degenerate by Corollary 
2.5.) rest point in D\. The following fact about Q± shows that our problem has fairly 
been solved. 



Proposition 4.4. With the above notations, Qi is a closed nonempty subset of Q 
with zero measure. 

Proof. We first show that Qi is closed and nonempty. Since D\ is compact and all 
rest points in the statements (A) and (B) of Theorem 4.1 are hyperbolic, both (A) 
and (B) occur in open subsets of Q. (In other words, our system is structurally stable 
in the nondegenerate case.) Since Q is connected, it cannot be the union of these 
two distinct open subsets. Thus Q\ is closed, nonempty and indeed large enough to 



separate two open subsets. We use Sard's theorem |L3[ to show that Q\ has zero 
measure. Notice that from (2-6), we can write i in terms of s if 7 + (£2 — £1 + A)s 7^ 
and from (2-7) we can write s in terms of i if (3\ + (A — e 2 )i 7^ 0. If Pi + (A — e 2 )i = 
then from i' = in (2-7), we have (e 2 — s\ — b — a) + (ei — £ 2 )i = 0. Thus 
e 2 — E\— 6 — a>0 and hence 7 + (e 2 — £1 + A)s > 0. It follows that either % can be 
written in terms of s from (2-6) or s in terms of i from (2-7). In the first case we 
have a root for the equation h\(s) = f3\ and in the latter case a root for h 2 (i) = 7 
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where: 

{% i 2 
hAs) = (b + Si - e 2 + a)- + (e 2 - X)i + (s 2 - £1) — , 
b s S is S 

h 2 {i) = — + (b + P)- + (A - £l )— - e 2 s, 

and r — 1 — s — i as before. Notice that in the first equation i has be written in 
terms of s and in the second equation s has been written in terms of i. Since at a 
degenerate rest point of the system (2-6), (2-7), the curves i' — and s' — are 
not transverse, it makes /3j a critical value of h\ or 7 a critical value of /i2- In other 
words, whenever the system (2-6), (2- 7) has a degenerate rest point in D\, either 0i 
or A belongs to the set of critical values which has zero measure by Sard's theorem. 
Therefore fli is contained in the union of two sets with zero measure. □ 



5 The Index Lemma 

A basic fact which has been used during the proof of our results is that the Poincare 
index of a piece-wise smooth Jordan curve on which the vector field is either tangent 
or inward is always 1. Here we provide the proof of this fact. The reader is referred 



to [|T^j for more details about the Poincare index. 

Let 7 : [0, 1] — ► M 2 be a piecewise smooth Jordan curve i.e. there exists a sequence 
= t < ti < . . . < t n = 1 such that 7 is smooth on for < % < n — 1. 

Suppose 7' is always nonzero on (ti,ti+i), moreover the left and right derivatives 
of 7 at ti exist and both are nonzero. With these assumptions we can define the 
external angle 0j G (— 7r,7r) at ti. Also the inward normal vector N(t) is defined for 
t ^ U and its right and left limit N(tf) exist at each U. For such a curve 7 we prove 
the following lemma. 

Lemma 5.1. Let U be a neighborhood of the image of 7 and X : U —> M 2 
be a smooth vector field which does not vanish on the image of 7 and satisfies 
N(t) ■ A( 7 (t)) > for t ^ U. Then I 7 (X) = 1. 
Proof. We define a new curve 71 : [0, n + 1] — > R 2 by 



7i (t) 




i + U <t <i + t i+1 , (0<z<n-l) 
i + ti+x < t < i + 1 + t i+ x- 



Notice that 71 moves like 7 but stops at 7(£j) for a unit of time, hence I-y(X) = 
I 71 (X). Now we use this unit of time to rotate N(t^) to arrive at N(tf). To do 
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this, we define a continuous function N\ : [0, n + 1] — > R 2 by 



N(t-i) i + ti<t<i + t i+1 , 
Rt0.(N(t-)) i + t i+1 <t<i + l + t i+1 , 



where R a is the rotation function with the angle a. Now both A"(7i(t)) and iVi(t) 
are continuous and do not vanish on [0, n + 1]. 

Claim: JVi(t) • Xfr^t)) > 0. 

Since iV(t) • X(i(t)) > for t ^ i*, by continuity we have A(£±) • X(7(^)) > 0. 
Moreover 0, e (— n, 7r), thus for any t £ [0, 1], there exists t' £ [0, 1] such that 



So i?^(A^(t7)).X(7(tj)) > 0. As a result of the above claim 

|A0(X( 7l (t))) - A0JVi(i)| < tt (5-1) 
where A6 is the total change of the angle on [0, n + 1]. Then we can write 
A0(X( 7l (t))) = A0(X( 7 (t))) = 2nI 1 (X). 

n—l n—l 

Ae(JVx(t)) = J] Aecjvxi^,^)) + Ae(Ar 1 | (i+ti+1 , i+1+ti+l) ) 

i=0 i=0 
n—l n 

i=0 i=l 

= 2tt. 

by Gauss-Bonnet theorem. Now from (5-1) we get \2ttI 1 (X) — 2n\ < it, hence 
\I 7 (X) — 1| < -. Since I 7 (X) is an integer, we get I-y(X) = 1. □ 
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